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Abstract 

Since it became evident that the CDM model for cosmic structure formation 
predicts smaller power on large scales than observed, many alternatives have 
been suggested. Among them, the existence of late decaying particle can cure 
it by delaying the beginning of the matter domination and increasing the 
horizon length at that time. We discuss the realization of this scenario and 
present the light neutrino and the light axino as possible examples of working 
particle physics model. We point out that the increased power at sub-galaxy 
scale predicted by this scenario could lead to rich sub-galaxy structures. 
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I. INTRODUCTION 



In the past decade, the cold dark matter (CDM) model has attracted a great attention 
as a possible theory of the cosmic structure formation. It is characterized by a spatially 
flat universe with ~90% of the mass density formed by cold dark matter and the scale 
invariant primordial fluctuation spectrum. It was supported by its successful features, the 
simplicity of its underlying assumptions and its link to the physics of the early universe. 
The basic idea that the large scale structure seen today evolved from very small primordial 
density fluctuations was strengthened by the recent detection of large scale anisotropies in 
the cosmic microwave background Q. The candidates for cold dark matter were provided 
by particle physics models. The scale invariant fluctuation spectrum was naturally obtained 
from inflation which is an indispensable ingredient of current cosmology. 

Recently, however, it became evident that this model is in conflict with observations 
0. The main difficulty was that the predicted power spectrum could not fit well with 
all the observational data simultaneously. The reason might be attributed to the relative 
smallness of the horizon length when the universe passes from relativistic matter domination 
to nonrelativistic matter domination which is the unique length scale characterizing the 
simple CDM model. Therefore, a more complicated model may be needed. 

There are two main variables in the structure formation theory, the content of matter 
in the universe and the shape of the primordial density fluctuations. A variation of two 
ingredients leads to a different prediction on the power spectrum. Actually several variants 
yield better agreements with observations than the cold dark matter model 0. The cold 
plus hot dark matter model attracted broad interest recently in this regard The warm 
dark matter model is another appealing possibility which draws new attention || . Altering 
the primordial fluctuation spectrum is another alternative . In addition, there is a way to 
maintain the matter content as cold dark matter by increasing the horizon length at the time 
of radiation matter equality by appropriate amount. Introduction of cosmological constant 
is one way to achieve this scenario J7|. The existence of late decaying particles which shift 
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the beginning of matter domination era is another possibility 0. At present, the large scale 
structure phenomenology does not favor any specific one among these variants. 

From the point of view of particle physicists, we have a similar situation because most 
dark matter candidates have no strong support for their existence from the particle physics 
phenomenology and variant models require unconfirmed new physics. Therefore, the present 
criterion in particle physics is the necessity and plausibility of the models accommodating 
these variants. In view of this criterion, the problem in the hot plus cold dark matter model is 
the one of explaining the comparable amount of hot and cold dark matter without attributing 
it to a mere chance. At first glance, the massive neutrino and the lightest supersymmetric 
particle as the hot and cold dark matter candidates seem to have no relation between them. 
Recently, some suggestions were made, connecting the amount of cold dark matter to that 
of hot dark matter which is usually assumed to be neutrinos || . For the cold plus hot dark 
matter models, one usually introduces additional parameters which need explanation. In a 
similar manner, introduction of late decaying particles faces the problem of explaining the 
simultaneous existence of cold dark matter with critical density and a late decaying particle 
with appropriate lifetime. It would face the similar problem of plausibility as the hot plus 
cold dark matter model does. At present, however, it is welcome to have a working model 
for it. This is our philosophy in this paper. We pursue for the possibility of late decaying 
particles dominating the energy density of the universe for some time after nucleosynthesis 
and decaying to extremely relativistic, weakly interacting particles. They shift the time of 
radiation matter equality, which leads to the elongation of the horizon length at that time 
and the suppression of the power on large scales compared to the value predicted in the 
simple CDM model. This idea was first pointed out in Ref. ||. Later, the 17 keV neutrino 
was used to implement the idea ||, but its existence is questioned now [[K]]. Recently, this 



idea has been rediscovered in the light axino decay [[□]] , through the naive expectation that 
the late decaying particles producing relativistic particles can mimick the cold plus hot dark 
matter scenario. 

In this paper, we emphasize the role of late decaying particles in cosmic structure for- 



mation and obtain the conditions on the relic amount and lifetime applicable to the generic 
case. Then the explicit models for late decaying particles are presented. The plan of the 
paper is as follows. In section II, we give a brief review on the structure formation theories 
and compare the several alternatives to the simple CDM model. In section III, the role of 
late decaying particles in structure formation is explained and the condition for better agree- 
ment with the observational data is presented. In section IV, we present the scenarios for 
the light neutrino and the light axino as explicit models for the cosmology of late decaying 
particles. Section V is a discussion and a conclusion. 



II. THE STRUCTURE FORMATION THEORIES AND DARK MATTER 

Structure formation and dark matter are closely related. Any theory for the formation of 
cosmic structure must specify two chief ingredients: the amount and nature of the material 
that fills the universe, and the properties of the seed density fluctuations from which galaxies 
and clusters developed. 

The amount of each component of materials in the universe is represented by the ratio 
of its mass density to the critical density: f2j = Pi/ p CT it- There are two widely accepted 
beliefs concerning the amount of matter [|l^]. One is that the total amount of matter is 
nearly critical (Q = Yli = 1)- The other is that the amount of baryons is less than 10% of 
the critical density (Qb ~ 0.1). The former follows from most inflationary scenarios. The 
strongest argument for the latter is the light element production in nucleosynthesis. The 
gap between Q and Qb can be filled by the existence of nonbaryonic dark matter or/and the 
cosmological constant. Dark matters are classified into hot, warm and cold dark matters by 
the characteristic scales they develop during the evolution of seed density fluctuations. 

At present, two types of origins for the seed density fluctuations are widely discussed: 
inflation and topological defects. In the inflationary scenario the quantum fluctuations of 
the inflaton field are changed into the density fluctuations, which is featured by the scale 
invariant fluctuation spectrum |0. The topological defects are created during the cosmic 
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phase transitions in some models with a spontaneously broken global symmetry. Cosmic 
strings and textures have been discussed in this context. 

The seed density fluctuations evolve under the act of gravity. For small fluctuations, the 
evolution is studied by the numerical integration of coupled perturbed-Einstein-Boltzmann 
equations. With the seed density fluctuation spectrum as initial condition, we obtain the 
evolved linear fluctuation spectrum. For large fluctuations, numerical simulation is used to 
show the nonlinear dynamics of clustering. Much work has been done with various dark 
matter contents and seed density fluctuation spectra. In the literature, one can find the 
numerical fit formulae for the evolved linear fluctuation spectra and the results of numerical 



simulations for many cases ||14|| . We can understand the qualitative features of these spectra 
by the characteristic scales inherent to the seed density fluctuations or developed by the 
matter content during the evolution. Important ones are the following. The Jeans length Aj 
separates the gravitationally stable and unstable modes. The free streaming length Xps is 
a scale below which fluctuations in a nearly collisionless component are damped due to free 
streaming. Another significant length scale is the horizon length at the time of the radiation 
matter equality Xeq- This scale appears because the growth of the density fluctuation of 
nonrelativistic matter within the horizon is suppressed during the radiation dominated era 
while begins as the matter domination era starts. 

Among the various models, the CDM model was most successful in the past decade. The 
CDM model assumes a flat (Q = 1) universe in which 5 ~ 10% of the critical density is 
provided by ordinary matter and the rest by the weakly interacting, yet unknown, cold dark 
matter. In addition, the seed fluctuations are assumed to be of the scale invariant form. 
The key assumptions of this model fit nicely with ideas on inflation. The evolved fluctuation 
spectrum is given by |L4] 

ir 12 _ ^ M 

1 kl (1 + ak + /3k 3 / 2 + 7 fc 2 ) 2 1 ' 

where A is a normalization constant and a = 1.7 I, (3 = 9.0 / 3//2 , 7 = 1.0 I 2 , and / = 
(Clh 2 )~ 1 6 1 / 2 M.pc. The scaling of I with (VLh 2 )~ l 6 1 / 2 simply reflects the horizon length at the 



time of radiation matter equality which is given by 

X EQ ~ 30(fi/i 2 )"^ 1/2 Mpc, (2) 

where 9 = p re i/1.68p 7 measures the present energy density of extremely relativistic particles 
and is 1 for photon and three massless neutrino flavors. It is independent of the nature 
of CDM and Xeq is a single physical length scale characterizing the above spectrum. The 



power spectrum of the model is compared to that extracted from the observational data |15 
in Fig. 1. It shows that with the standard choice flh = 0.5 the model does not fit the data. 
With the normalization fixed at COBE data, we can say that the standard CDM model 
predicts more power at small scales than observed. From the figure, we intuitively reason 
that the problem in the CDM model is that it has a relatively small A^q. 

Several cures for the excessive power of the simple CDM model at small scales have been 
considered. We put them into four categories: (i) assuming large biasing or antibiasing, (ii) 
the initial fluctuation spectrum with less power at small scales, (iii) the mixed dark matter 
contents, and (iv) increasing Xeq- Biasing is the modulation of the galaxy distribution by 
processes associated with galaxy formation. Though a complicated version in this category 
met some successes, we don't know its detailed physical mechanism. The initial fluctuation 
spectrum with less power at small scales can be obtained in fine tuned inflationary models 
H or in cosmic string and texture models. But they require peculiar models with fine tuning 
of parameters. The mixed dark matter content will change the evolution of the initial power 
spectrum. Cold plus hot dark matter model was shown to be successful in this regard. 
With the massive neutrino as hot dark matter, the best fit was obtained for ^cdm — 0.6, 
VL v = 0.3 and fie = 0.1 |16| . While this model met many cosmological successes, it also 
has a plausibility problem: how to explain the comparable amount of hot and cold dark 
matter without attributing it to a mere chance. It might lead to fine tuning again. Warm 
dark matter model may give a promising result. In increasing Xeq, the simplest idea is 
to lower Q or h to make Qh ~ 0.2 — 0.3. (See Fig. 1. We allow small biasing.) But this 
conflicts with other observational data JL^]. Introduction of small cosmological constant 



allows Ocdm = 0.2 in the flat universe, resulting in the increase of Xeq 0. But this requires 
severe fine tuning of the cosmological constant, which is not solved yet in particle physics. 
Late decaying particles can lead to the retardation of matter domination with the desired 
increase in Xeq- We will discuss this scenario in detail in the following section. 

III. THE ROLE OF LATE DECAYING RELIC PARTICLES IN STRUCTURE 



In the CDM model, the power on small scales can be reduced by delaying the beginning of 
matter domination era. Since the density fluctuations of cold dark matter grow significantly 
during the matter domination era, the growing time for small scales is reduced as much as 
one delays the time for the small scales to enter the horizon before the radiation matter 
equality. The Qh = 0.2 CDM model is very successful in this way. But the problem is that 
one wants to get the same success while maintaining Q — 1. One way to achieve this is by 
postulating a late decaying relic particle which dominates the energy density of the universe 
for some time and decays into extremely relativistic particles. By producing relativistic 
particles, the universe enters into a radiation dominated era around the time of its decay. 

Consider a late decaying relic particle X with the mass m, the lifetime r ^> 1 sec (the 
time of neutrino decoupling) and the ratio of the number density to entropy density Y 
which is kept constant when the processes creating or destroying X are frozen. The value 
of Y varies from ~ 10~ 2 to 10 -9 . The number density of X at temperature T is given 
by Tlx CO = Ys(T), where s(T) = ^r-g* s (T)T 3 is the entropy density. For T C m, 1 is 
nonrelativistic and the energy density is given by px{T) = mnx{T). 

Let us consider the relevant range of the values of mY and r. There are two obvious 
bounds on mY and r. The first concerns nucleosynthesis. At the time of neutrino decoupling, 
one customarily requires that the energy density of X should be less than that of one neutrino 
species. This imposes an upper bound on mY which applies to the case that m > 1 MeV, 
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The second comes from the condition that the energy density of the decay products of X 
should not overdose the universe: 

f-K^ 2 <2xl0 6 /, 3 . (4) 
Vsec/ VMeVy v ; 

If X can decay into photons or ordinary particles, a severe restriction can be imposed by the 

present observation of microwave background radiation [[L7|j. Therefore, we consider cases in 



which the branching ratio of X into photons is negligible and decay products of X are very 
weakly interacting so as not to disturb the spectra of observable particles. 

In this scenario, there exist two eras of matter domination. The first era begins when the 
energy density of the late decaying particles dominates the energy density of the universe. 
The first matter dominated era ends when the late decaying particle decays into extremely 
relativistic particles. The second era begins when the energy density of cold dark matter 
dominates over the energy density of the decay products. The evolution of the energy 
densities is shown in Fig. 2. Four interesting scale factors appear in this scenario: Reqi 
(the scale factor at which the energy density of X's and the energy density of photons and 
neutrinos become equal), Re (the scale factor at which X decays), Req (the scale factor 
at which the energy density of cold dark matter and the energy density of photons and 
neutrinos become equal), and Reqi (the scale factor at which energy densities of cold dark 
matter and extremely relativistic decay products become equal). Req is the scale factor at 
the radiation-matter equality epoch in the simple CDM model. In our scenario, we make 
Req2 > Req- 

The temperature and the time at which the first matter domination begins are given by 

T EQ i = 1.55 mY, (5) 
/MeV\ 2 

t EQ i = 0.55 [—^J sec, (6) 

if the lifetime r is larger than t E Qi- Actually, the conditions for which X-domination really 
occurs are r > tEQi and t E Qi < tEQi which are converted to 



t \ / mY \ 



2 



sec/ VMeV 



J >0.55, (7) 
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To find out the scale factor Req2, we use the rough approximation of simultaneous decay 
and sudden change of radiation and matter domination. The condition for Req and Req2 
read 

\^eqJ \ReqJ 



( Ro V 
Po -5 — Peqi 



Reqi\ 3 ( Rd \ ( Reqi\ 



(10) 



Rd J \ Req2 J \ Req2 J 

where Peqi is the energy densities of X at the first radiation matter equality and po is the 
present matter energy density. From these, it follows that the second matter domination in 
our scenario occurs later than the well-known matter domination of the simple CDM model 
by a factor 



Req2 Rd 



2/3 



1~ — +1, (11) 



Req Reqi \tEQi, 

where the last equality comes from the relation R oc t 2//3 during the matter dominated era. 

The evolution of the fluctuation spectrum is obtained by the linear perturbation theory 
as done in Ref. || for the 17 keV neutrino case. But here, we present a rough behavior and a 
simple estimate. The fluctuation spectrum is characterized by two length scales, associated 
with the horizon lengths at two instances of radiation matter equality, 

9 /MeV\ , . , 

X EQ1 * 8 x 10- 2 (—J kpc, (12) 



X EQ 2 ^ S0(Qh 



2\-l 



1 f t \ ( rnY \ 
VMeV/ 



0.55 Vsec 



2^, 2/3 



+ 1 



1/2 



Mpc. (13) 



The existence of the scale Xeqi distinguishes this scenario from the standard CDM model. 
Objects on this and smaller scales would collapse at high red shifts and more structures exist 
on these scales than the standard CDM model. We will discuss about it later. 

Xeq2 is larger than Xeq when the X-domination condition, eq. (|7|) holds. In the previous 
section, we mentioned that the Qh = 0.2 CDM model have a good fit with possible small 
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biasing taken into account. Lowering VLh has an effect of making Xeq large. But this 
model is in conflict with the inflationary paradigm = 1 and the observed Hubble constant 
ft, ~ 0.4 — 1. In the late decaying particle scenario, we have the same effect of making Xeq 
large but still can keep Q — 1. The value flh = 0.2 gives a rough estimation of the desired 
value for \eq2- ^EQ2/h~ 1 Mpc ~ 150. It can be rewritten as a condition on r and mY for 
VL = 1, 



This together with the eqs. @ and (H) forms a main condition of the late decaying particle 
scenario with £1 = 1. 



The long lifetime and the small radiative branching ratio of a late decaying particle make 
it difficult to construct particle physics models for it. But it is certainly not impossible. It 
is possible to use symmetry and/or coupling suppressed by a large mass scale to obtain the 
long lifetime and to suppress the radiative decay. In this section, we present two models: 
the massive neutrino and the light axino. 



Cosmological implications of the massive neutrino vary depending on the size of mass 
and lifetime. At present, the stable ~ 30 eV mass neutrino is a good candidate for hot dark 
matter. Here we consider a different region in mass and lifetime. The relic abundance of a 
light neutrino species is Y = 3.9 x 10~ 2 . The late decaying particle scenario works for the 
neutrino mass in the range 100 eV ^ m v < 1 MeV. Hereafter, we use m v = 17 keV as the 
representative value and adopt the notation run = m u /17keV . Then mY = 0.6Qmn keV 
and the eq. (|T4]) demands the lifetime to be 




(14) 



IV. MODELS OF A LATE DECAYING PARTICLE 



A. A Late Decaying Neutrino 



r ~ 0.5m 17 2 yr, 



(15) 
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for the standard value h = 0.5. (Actually, the use of the formulae of the previous sec- 
tion needs a slight modification due to the change in the number of relativistic neutrino 
species. But we still use them because the estimations are not affected very much by this 
modification.) The length scale of the first radiation matter equality is 

\ EQ1 ~ 120m ^ kpc. (16) 

The well-known way to obtain the neutrino mass in the above mentioned range is to use 
the seesaw mechanism. The long lifetime required in the eq. ( |l5l) can be obtained in several 
models fl8 |. They considered schemes in which the massive neutrino decays into exotic 



particles — majoron, familon or techniphoton. In either case, we have to make a judicious 
choice of unknown parameters to give such a long lifetime. 



B. A Light Axino and A Lighter Gravitino 

The importance of the cosmological effects of axinos and gravitinos arises from the weak- 
nesses of their interactions due to the suppression factors of the axion decay constant for 
the axino and the Planck mass for the gravitino. Here we focus on the scenario in which the 
gravitino is the lightest supersymmetric particle (LSP) and the axino is the second lightest 
super symmetric particle (SLSP). 

For this, we need the simultaneous implementation of local supersymmetry and the axion 
solution of the strong CP problem. Supersymmetry has been introduced to solve the gauge 
hierarchy problem. Realistic models incorporating low energy supersymmetry consider local 
supersymmetry (supergravity), spontaneously broken by the hidden sector [fL9| . Most of 
them have an R symmetry to suppress the unwanted large B and L violations. The R 
symmetry dictates the existence of a stable particle commonly called the LSP which is a 
good candidate for dark matter. The LSP is model dependent and a popular candidate 
is neutralino, a linear combination of neutral gauginos and neutral higgsinos. If we adopt 
local supersymmetry, the gravitino is another natural candidate for the LSP. Peccei-Quinn 
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symmetry is needed for a solution of the strong CP problem, and no idea seems to be 



compelling in this regard ||20|. It predicts the existence of the axion [21 1 and the coherent 



classical axion oscillation created at T ~ Aqcd during the expansion of the universe is a 



good candidate for dark matter [22]. The simultaneous implementation of these two ideas 



can affect the status of cold dark matter candidates P5| . The axino (the superpartner of 



axion) can and generally do destabilize the lightest ordinary supersymmetric particle (LOSP) 
which is a popular dark matter candidate. Then the viable candidates for dark matter are 
gravitinos, axinos and axions. For the gravitino and the axino to be LSP and SLSP, we need 
gravitino and axino lighter than the LOSP. 

Currently favored gravitino mass is identified with the electroweak scale of 10 2 — 10 3 
GeV and the gravitino with such mass cannot be the LSP because there is usually a lighter 
supersymmetric particle to which the gravitino can decay. Recently, however, the old idea of 
supercolor at multi TeV region has been reinvestigated without phenomenological difficulties 
||24|| . In this case the gravitino mass falls in the 0.4 meV — 0.4 keV region for A superco i or = 
1 — 10 3 TeV. Therefore, the phenomenologically favored gravitino mass is still open in these 
two regions and it is worthwhile to consider the cosmological consequence of this light 
gravitino scenario. 

In the case that R symmetry is exact and the axino is lighter than the LOSP, the LSP 
lighter than the axino must exist for the axino to decay. The presently known candidate is 
the gravitino. The possibility and cosmological constraints of the gravitino as the LSP were 
previously considered [p5|| . Here, we summarize the results. Gravitinos typically decouple 
just after the Planck time. From this one obtains a gravitino mass bound 7723/2 < 2 keV 
in the standard big bang cosmology. But the inflation after the gravitino decoupling might 
have washed out relic gravitinos almost completely and we have no such bound. Gravitinos 
can be regenerated through the reheating process subsequent to inflation. The requirement 
that the regenerated gravitinos should not overdose the universe gives a constraint between 
the reheating temperature and the gravitino mass PS|. If we assume that the reheating 



temperature is lager than the axino decoupling temperature 10 11 GeV, a region of gravitino 
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mass from a few keV to a few x 10 GeV is excluded. There is also a bound in the extremely 
light gravitino. For the very light gravitino whose mass m 3 / 2 is small compared to the typical 
interaction energy, the helicity ±| component of the gravitino dominate the interaction and 
the strength of the interaction becomes larger as m 3 / 2 becomes smaller. Therefore, the 
extremely light gravitino (m 3 / 2 <C 1 keV) interacts strongly and decouples well below the 
weak scale. The standard nucleosynthesis scenario constrains the number of the species 
of the neutrino-like particles to be smaller than 3.3 at the time of nucleosynthesis. Then, 
gravitinos must decouple before T ~ 200 MeV to get a sufficient dilution factor. This yields 
a gravitino mass bound 

m ^ io - 4 (io5y ev < i7 > 

where mj is the slepton mass. 

The axino mass is predicted to be in a region between 10 GeV and keV in supergravity 
models with supersymmetry broken in the hidden sector. Usually one adopts the minimal 
kinetic energy term. In general, there exist other loop contributions and the kinetic energy 
term can be of nonminimal form. Furthermore, we can adopt the supersymmetry breaking 
by supercolor type interactions. Therefore, the axino mass heavier than keV is a theoretical 
possibility. 

Now, we calculate the axino lifetime. The dominant coupling of axion to gravitino is 

£aaG = ^,l U d u Z^~a L + h.C. (18) 

where z = (s + ia)/\/2 in terms of saxino s and axion a. There exists another interaction 
term proportional to m 3 / 2 which is negligible compared to the above interaction. In the case 
that the axino is heavier than the gravitino, the axino decays into the gravitino and the 
axion. In the limit of vanishing axion mass, the axino lifetime is 

= 9Mt yi /2 = L2 x i()12 ,MeVy ^ 
ra~ V m~ a ) V eV / 

Axinos are kept in thermal equilibrium at high temperature through the process qq <-> ag 

with an intermediate gluon. From the interaction lagrangian 
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C = £k ( aG ^ da " V + sG % Qa " u + ~^,u~g a G a ^ + •••), (20) 

the reaction rate can be calculated to give V ~ ct 3 T 3 /16nF 2 for T < F a . On the other hand, 
the expansion rate is H = l.G6g¥ 2 T 2 /M p i where g* = 915/4 for the particle content in the 
minimal supersymmetric standard model. Therefore, the decoupling temperature is 

" n (lo4v) 2 © 3GeV - < 21 > 

The amount of relic axinos depend on the decoupling temperature. (There are also 
axinos from decays of relic LOSPs. But the number density of relic LOSPs is usually much 
smaller than that of hot thermal relic and hence axinos from decays of relic LOSPs can be 
neglected.) The axino decoupling temperature T a is so high that it can be comparable to 
or larger than the reheating temperature Tr. For simplicity, we consider two extreme cases 
which differ from each other significantly. 

(i) Tr ^> To,: The relic axinos are hot thermal relic so that Y ~ 1.8 x 10~ 3 . The desired 
value of \eq2 is obtained for 



within the mass range 

2 keV < m- a < 60 MeV, (23) 

which come from the eqs. (|3]), (§) and @. This requires the MeV mass axino and the 
very light, sub-eV mass gravitino. The cosmologically allowed mass ranges and the region 
corresponding to the eq. fl22|) are shown in Fig. 3. The smaller characteristic length scale is 
given by 

Xeqi = 44 (^) kpc. (25) 
With the MeV mass axino, it is about the globular cluster scale. 
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(ii) Tr <C Ta~. The hot thermal relic axinos are diluted away by inflation and relic axinos 
are those regenerated through the nonequilibrium processes, which yields Y ~ 7x 10~ 6 F~ 2 Tr 
where F a = F a /10 12 GeV and Tr = Tjj/10 6 GeV. Now the corresponding equations of the 
eqs. (p2j)-(|25|) contain F a , Tr as well as ma, 777,3/2 and are given by 



^(^) 3 (^) 2 --«(- 2 - 1 



3/2 



(26) 



5 x 10- 4 < F~ 2 Tr ( j < 15, (27) 



^(^) 3 (^) 2 <-^ ! 



(28) 



and 



A^^llF^ 1 ^) 2 kpc. (29) 
\ ilia / 

This case opens a possibility of the heavier axino with GeV mass depending on the size of 
the reheating temperature. 

V. DISCUSSION AND CONCLUSION 

In our scenario of late decaying particles, there exists a brief era of matter domination 
between Xeqi and X^. During this epoch, there exists a possibility that some scales can 
grow and form structures. Because there exists a characteristic scale of globular clusters, 
it is worthwhile to see if the structure formation is triggered during this brief period and 
to check whether it coincides with the scale of globular clusters. If the conclusion is af- 
firmative, then the scenario of the late decaying particle provides an explanation for the 
size of globular clusters. Furthermore, such a large power at scale smaller than the galaxy 
scale is advantageous on explaining the Lyman-line-absorbing clouds seen in quasar spectra 



1 27]. The problem is that it leads to large anisotropies in microwave background radiation 
at small angles. But it is known that this early structure formation, e.g. the early pop- 
ulation of quasi-stellar objects, could reionize the universe shortly after the normal epoch 
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of photon-baryon decoupling at z ~ 1000, thereby reducing small angle anisotropies in the 
microwave background radiation. 

There are a few merits of the axino-gravitino scenario. First, it has the classical axion 
oscillation as a natural candidate for cold dark matter. Second, the interactions of the axino 
and the gravitino are well known so that the unknown parameters are mainly their masses 
and no unknown coupling is involved. Third, the explanation of the globular cluster scale 
or other sub-galaxy structures may fix their masses. 

In summary, we considered the late decaying particle scenario in the Q = 1 CDM model 
as a viable alternative to the simple CDM model which failed in fitting all the observational 
data simultaneously. The late decaying particles delay the beginning of the matter domi- 
nation, thereby reducing the excessive power at small scales which was the problem of the 
simple CDM model. This scenario has another peak in the power spectrum developed during 
the earlier matter domination. Its length scale is smaller than the galaxy scale, thereby hav- 
ing implications for sub-galaxy structures, which need further investigations. The particle 
physics model realizing this scenario is difficult to construct due to the required long lifetime 
and small radiative branching ratio of a late decaying particle, but certainly not impossible. 
We presented two examples. The keV mass neutrino with lifetime of order 1 yr, which has 
been considered by many authors, is a possible candidate of the late decaying particle. The 
axino-gravitino scenario requires the MeV mass axino and the sub-eV mass gravitino (the 
sufficiently low reheating temperature can raise the required masses). It possesses all the 
necessary ingredients of the late decaying particle scenario with the simultaneous implemen- 
tation of local supersymmetry and Peccei-Quinn symmetry. 
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FIGURES 

FIG. 1. Comparison of the CDM model power spectrum with the observational data. The solid 
line corresponds to the £lh = 0.3 CDM model, while the dashed and dotted lines correspond to the 
Qh = 0.5 and 1 CDM models respectively. 

FIG. 2. The evolution of energy densities of radiation (7, u), late decaying particles (X), their 
decay products and the cold dark matter. 

FIG. 3. Cosmologically allowed region in the axino mass - gravitino mass plane. Dotted regions 
are excluded by the closure bound and the nucleosynthesis bound (the equation numbers represent 
the relevant equations in the text). The heavy dotted line corresponds to the eq. (p2|) for h = 0.5. 



The white island below the dashed line is the mass range allowed by the eqs. (23) and fl24|). 



22 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405385v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405385v2 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405385v2 



